arXiv:astro-ph/0503065vl 3 Mar 2005 


February 2, 2008 16:44 Proceedings Trim Size: 9in x 6in 


otasi4up 


astro-ph/0503065 
UMN-TH-2346/05 
FTPI-MINN-05/05 


March 2005 


TASI LECTURES ON ASTROPARTICLE PHYSICS 


KEITH A. OLIVE* 

William I. Fine Theoretical Physics Institute, School of Physics and Astronomy, 
University of Minnesota, Minneapolis, MN 55455 USA 
E-mail: olive@umn.edu 


Selected topics in Astroparticle Physics including the CMB, dark matter, BBN, 
and the variations of fundamental couplings are discussed. 


1. Introduction 

The background for all of the topics to be discussed in these lectures is 
the Big bang model. The observed homogeneity and isotropy enable us 
to describe the overall geometry and evolution of the Universe in terms of 
two cosmological parameters accounting for the spatial curvature and the 
overall expansion (or contraction) of the Universe. These two quantities 
appear in the most general expression for a space-time metric which has 
a (3D) maximally symmetric subspace of a 4D space-time, known as the 
Robertson-Walker metric: 



( 1 ) 


where R(t) is the cosmological scale factor and k is the curvature constant. 
By rescaling the radial coordinate, we can choose k to take only the dis¬ 
crete values +1, —1, or 0 corresponding to closed, open, or spatially flat 
geometries. 

The cosmological equations of motion are derived from Einstein’s equa¬ 
tions 


( 2 ) 




*This work was supported in part by DOE grant DE-FG02-94ER40823 at Minnesota. 


1 












February 2, 2008 16:44 Proceedings Trim Size: 9in x 6in 


otasi4up 


2 


where A is the cosmological constant. It is common to assume that the 
matter content of the Universe is a perfect fluid, for which 


= -pg^u + {p + p) u^Uv 


(3) 


where g is the space-time metric described by (1), p is the isotropic 
pressure, p is the energy density and u = (1, 0, 0,0) is the velocity vector for 
the isotropic fluid in co-moving coordinates. With the perfect fluid source, 
Einstein’s equations lead to the Friedmann-Lemaitre equations 


and 
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where Hit) is the Hubble parameter. Energy conservation via T// = 0, 
leads to a third useful equation [which can also be derived from Eqs. (4) 
and (5)] 


p = -3H(p + p) (6) 

The Friedmann equation can be rewritten as 

(n - i )h 2 = A (7 ) 

so that k = 0,+1,—1 corresponds to ft = l,ft > 1 and SI < 1. However, 
the value of SI appearing in Eq. (7) represents the sum SI = ft m + St a of 
contributions from the matter density (St m ) and the cosmological constant 
(St A = A/377 2 ). 


2. The CMB 

There has been a great deal of progress in the last several years concerning 
the determination of both St m and St A . Cosmic Microwave Background 
(CMB) anisotropy experiments have been able to determine the curvature 
(i.e. the sum of St m and St A ) to within a few percent, while observations 
of type la supernovae at high redshift provide information on a (nearly) 
orthogonal combination of the two density parameters. 

The CMB is of course deeply rooted in the development and verification 
of the big bang model and big bang nucleosynthesis (BBN) 1 . Indeed, it was 
the formulation of BBN that led to the prediction of the microwave back¬ 
ground. The argument is rather simple. BBN requires temperatures greater 
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than 100 keV, which according to the standard model time-temperature re- 
lation, fs^Mev = 2A/y/N, where N is the number of relativistic degrees of 
freedom at temperature T, and corresponds to timescales less than about 
200 s. The typical cross section for the first link in the nucleosynthetic 
chain is 

av(p + n —> D + 7)~5x 10~ 20 cm 3 /s (8) 

This implies that it was necessary to achieve a density 

n ~-~ 10 17 cm -3 (9) 

avt 

for nucleosynthesis to begin. The density in baryons today is known ap¬ 
proximately from the density of visible matter to be n Bo ~ 10~ 7 cm -3 and 
since we know that that the density n scales as R~ a ~ T 3 , the temperature 
today must be 

T 0 = (n Bo /n) 1/3 T B bn ~ 10K (10) 

thus linking two of the most important tests of the big bang theory. 

An enormous amount of cosmological information is encoded in the 
angular expansion of the CMB temperature 

0 )* ( 11 ) 

The monopole term characterizes the mean background temperature of 
T 7 = 2.725 ± 0.001 K as determined by COBE 2 , whereas the dipole term 
can be associated with the Doppler shift produced by our peculiar motion 
with respect to the CMB. In contrast, the higher order multipoles, are 
directly related to energy density perturbations in the early Universe. When 
compared with theoretical models, the higher order anisotropies can be used 
to constrain several key cosmological parameters. In the context of simple 
adiabatic cold dark matter (CDM) models, there are nine of these: the 
cold dark matter density, fl x /i 2 ; the baryon density, fl B h 2 ; the curvature - 
characterized by fitotai! the hubble parameter, h; the optical depth, r; the 
spectral indices of scalar and tensor perturbations, n s and np, the ratio of 
tensor to scalar perturbations, r; and the overall amplitude of fluctuations, 

Q. 

Microwave background anisotropy measurements have made tremen¬ 
dous advances in the last few years. The power spectrum 3 ’ 4,5,6 ’ 7 ’ 8 ’ 9 ’ 10 has 
been measured relatively accurately out to multipole moments correspond¬ 
ing to £ ~ 2000. A compilation of recent data is shown in Fig. 1 11 , where 
the power in at each £ is given by (2£ + 1)Q/(47t), and Ci =< |a^ m | 2 >. 



February 2, 2008 16:44 Proceedings Trim Size: 9in x 6in 


otasi4up 


4 


Angular Scale 

1° 10' 5' 



Multipole l 


Figure 1. The power in the microwave background anisotropy spectrum as measured 
by WMAP 9 , CBI 6 , and ACBAR 10 . Taken from 11 . 


As indicated above, the details of this spectrum enable one to make 
accurate predictions of a large number of fundamental cosmological param¬ 
eters. The results of the WMAP data (with other information concerning 
the power spectrum) is shown in Table 1. For details see ref. 9 . 



WMAP alone 

WMAPext ± 2dFGRS 

WMAPext ± 2dFGRS ± Lyman a 


power-law 

power-law 

running 

n m h 2 

0.14 ±0.02 

0.134 ±0.006 

0 1 QC+U.0U8 

u.ioo_o 009 

n B h 2 

0.024 ±0.001 

0.023 ±0.001 

0.0224 ±0.0009 

h 

0.72 ±0.05 

0.73 ±0.03 

0 71 +0.04 
u -‘ -*--0.03 

n s 

0.99 ±0.04 

0.97 ±0.03 

0.93 ±0.03 

T 

0.i66±S;g;S 

0.148±g;g;? 

0.17 ±0.06 


Of particular interest to us here is the CMB determination of the total 
density, fltot, as well as the matter density fi m . There is strong evidence 
that the Universe is flat or very close to it. The best constraint on fl to tai is 
1.02 ±0.02. Furthermore, the matter density is significantly larger than the 
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baryon density implying the existence of cold dark matter and the baryon 
density, as we will see below, is consistent with the BBN production of D/H 
and its abundance in quasar absorption systems. The apparent discrepancy 
between the CMB value of fitot and fi m , though not conclusive on its own, is 
a sign that a contribution from the vacuum energy density or cosmological 
constant, is also required. The preferred region in the fi m — plane is 
shown in Fig. 2 under four different assumptions 9 . 



Figure 2. Two-dimensional confidence regions in the (Qm,^a) plane 9 . 


The presence or absence of a cosmological constant is a long standing 
problem in cosmology. We know that the cosmological term is at most a 
factor of a few times larger than the current mass density. Thus from Eq. 
(4), we see that the dimensionless combination, G^A < 10 -121 . Never¬ 
theless, even a small non-zero value for A could greatly affect the future 
history of the Universe: allowing open Universes to recollapse (if A < 0), 
or closed Universes to expand forever (if A > 0 and sufficiently large). 

When the SN la results 12 are included (see the last panel of Fig. 2) we 
are led to a seemingly conclusive picture. The Universe is nearly flat with 
fitot — 1- However the density in matter makes up only 23% of this total, 
with the remainder in a cosmological constant or some other form of dark 
energy. 
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3. Dark Matter 

3.1. Observational Evidence 

Direct observational evidence for dark matter is found from a variety of 
sources. On the scale of galactic halos, the observed flatness of the rotation 
curves of spiral galaxies is a clear indicator for dark matter. There is also 
evidence for dark matter in elliptical galaxies, as well as clusters of galaxies 
coming from the X-ray observations of these objects. Also, direct evidence 
has been obtained through the study of gravitational lenses. 

For example, assuming that galaxies are in virial equilibrium, one ex¬ 
pects that one can relate the mass at a given distance r, from the center of 
a galaxy to its rotational velocity by 

M{r) oc v 2 t/Gn (12) 

The rotational velocity, v. is measured 13,14 by observing 21 cm emission 
lines in HI regions (neutral hydrogen) beyond the point where most of 
the light in the galaxy ceases. A subset of a compilation 15 of nearly 1000 
rotation curves of spiral galaxies is shown in Fig. 3. The subset shown is 
restricted to a narrow range in brightness, but is characteristic for a wide 
range of spiral galaxies. Shown is the rotational velocity as a function of 
r in units of the optical radius. If the bulk of the mass is associated with 
light, then beyond the point where most of the light stops, M would be 
constant and v 2 oc 1/r. This is not the case, as the rotation curves appear 
to be flat, i.e., v ~ constant outside the core of the galaxy. This implies 
that M oc r beyond the point where the light stops. This is one of the 
strongest pieces of evidence for the existence of dark matter on galactic 
scales. Velocity measurements indicate dark matter in elliptical galaxies as 
well 16 . For a more complete discussion see 17 . 


3.2. Theory 

Theoretically, there is no lack of support for the dark matter hypothesis. 
The standard big bang model including inflation almost requires f2 t ot = 
1 18 . This can be seen from the following simple solution to the curvature 
problem. The unfortunate fact that at present we do not even know whether 
n is larger or smaller than one, indicates that we do not know the sign of 
the curvature term further implying that it is subdominant in Eq. (4) 

k 87 tG 

D2 ^ O P 


(13) 
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<lft> = -21.2 



R/Ropt 


Figure 3. Synthetic rotation curve 15 for galaxies with (M) = —21.2. The dotted curve 
shows the disk contribution, whereas the dashed curve shows the halo contribution. 


In an adiabatically expanding Universe, R ~ T 1 where T is the tempera¬ 
ture of the thermal photon background. Therefore the quantity 


k = 


< ^5 < 2 x 10~ 58 


R 2 T 2 3 T 2 


(14) 


is dimensionless and constant in the standard model. This is known as 
the curvature problem and can be resolved by a period of inflation. Before 
inflation, let us write R = Ri, T = Tt and R ~ T -1 . During inflation, 
R ~ T -1 ~ e Ht , where H is constant. After inflation, R = Rf Ri but 
T = Tf = Tr < Ti where Tr is the temperature to which the Universe 
reheats. Thus R ^ T and k — > 0 is not constant. But from Eqs. (7) and 
(14) if k —> 0 then 0 —> 1, and since typical inflationary models contain 
much more expansion than is necessary, Q becomes exponentially close to 
one. 

The existence of non-baryonic dark matter can be immediately inferred 
from the determination of the cosmological parameters through the mi¬ 
crowave background anisotropy as described above. If ~ 0.13 and 

n B h 2 ~ 0 . 02 , then the difference must be dark matter which contributes 
to the total density Odm ^ 2 — 0.11. In addition, because the amplitude of 
fluctuations is relatively small, dark matter is necessary to have sufficient 
time to grow primordial perturbations into galaxies (for a more complete 
discussion see 17 ). 
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3.3. Candidates 
3.3.1. Baryons 

Accepting the dark matter hypothesis, the first choice for a candidate 
should be something we know to exist, baryons. Though baryonic dark 
matter can not be the whole story if fi m > 0 . 1 , the identity of the dark 
matter in galactic halos, which appear to contribute at the level of f2 ~ 0.05, 
remains an important question needing to be resolved. A baryon density of 
this magnitude is not excluded by nucleosynthesis. Indeed we know some 
of the baryons are dark since fl < 0.01 in the disk of the galaxy. 

It is interesting to note that until recently, there seemed to be some dif¬ 
ficulty in reconciling the baryon budget of the Universe. By counting the 
visible contribution to fi in stellar populations and the X-ray producing hot 
gas, Persic and Salucci 19 found only fl v i s ~ 0.003. A subsequent account¬ 
ing by Fukugita, Hogan and Peebles 20 found slightly more (fl ~ 0.02) by 
including the contribution from plasmas in groups and clusters. At high 
redshift on the other hand, all of the baryons can be accounted for. The 
observed opacity of the Ly a forest in QSO absorption spectra requires a 
large baryon density consistent with the determinations by the CMB and 
BBN 21 . 

In galactic halos, however, it is quite difficult to hide large amounts of 
baryonic matter. Sites for halo baryons that have been discussed include 
Hydrogen (frozen, cold or hot gas), low mass stars/Jupiters, remnants of 
massive stars such as white dwarfs, neutron stars or black holes. In almost 
every case, a serious theoretical or observational problem is encountered 22 . 


3.3.2. Neutrinos 

Light neutrinos (in < 30eU) are a long-time standard when it comes to 
non-baryonic dark matter 23 . Light neutrinos are, however, ruled out as a 
dominant form of dark matter because they produce too much large scale 
structure 24 . Because the smallest non-linear structures have mass scale 
Mj « 3x lO 18 Af 0 /?Ti 2 (eV) and the typical galactic mass scale is ~ 10 12 -M Q , 
galaxies must fragment out of the larger pancake-like objects. The problem 
with such a scenario is that galaxies form late 25,26 (z < 1 ) whereas quasars 
and galaxies are seen out to redshifts z > 6. 

The neutrino decoupling scale of 0(1 ) MeV has an important conse¬ 
quence on the final relic density of massive neutrinos. Neutrinos more 
massive than 1 MeV will begin to annihilate prior to decoupling, and while 
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Figure 4. Summary plot 28 of the relic density of Dirac neutrinos (solid) including a 
possible neutrino asymmetry of ??„ = 5 X 10“ 11 (dotted). 

in equilibrium, their number density will become exponentially suppressed. 
Lighter neutrinos decouple as radiation on the other hand, and hence do 
not experience the suppression due to annihilation. Therefore, the calcu¬ 
lations of the number density of light (to„ < 1 MeV) and heavy (m„ > 1 
MeV) neutrinos differ substantially. 

The energy of density of light neutrinos with m v < 1 MeV can be 
expressed at late times as = m v Y u n 1 where Y v = n^/n 7 is the number 
density of p’s relative to the density of photons, which today is 411 photons 
per cm 3 . It is easy to show that in an adiabatically expanding universe 
Y v = 3/11. This suppression is a result of the e + e _ annihilation which 
occurs after neutrino decoupling and heats the photon bath relative to the 
neutrinos. Imposing the constraint Q„h 2 < 0.13, translates into a strong 
constraint (upper bound) on Majorana neutrino masses 27 : 

mtot = T, < 12eV. (15) 

V 

where the sum runs over neutrino mass eigenstates. The limit for Dirac 
neutrinos depends on the interactions of the right-handed states. The limit 
(15) and the corresponding initial rise in £l u h 2 as a function of m v is dis¬ 
played in the Figure 4. 

Combining the rapidly improving data on key cosmological parameters 
with the better statistics from large redshift surveys has made it possible 
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to go a step forward along this path. It is now possible to set stringent 
limits on the light neutrino mass density fand hence on the neutrino 
mass based on the power spectrum of the Ly a forest 29 , m to t < 5.5 eV, 
and the limit is even stronger if the total matter density, fl m is less than 
0.5. Adding additional observation constraints from the CMB and galaxy 
clusters drops this limit 30 to 4.2 eV. This limit has recently been improved 
by the 2dF Galaxy redshift 31 survey by comparing the derived power spec¬ 
trum of fluctuations with structure formation models. Focussing on the 
the presently favoured ACDM model, the neutrino mass bound becomes 
mtot < 1-8 eV for fi m < 0.5. When even more constraints such as HST 
Key project data, supernovae type la data, and BBN are included 32 the 
limit can be pushed to mtot < 0.9 eV. With WMAP data, an upper limit 
of m to t < 0.7 eV has been derived 9 . 

The calculation of the relic density for neutrinos more massive than ~ 1 
MeV, is substantially more involved. The relic density is now determined 
by the freeze-out of neutrino annihilations which occur at T < m u , after 
annihilations have begun to seriously reduce their number density 33 . For 
particles which annihilate through approximate weak scale interactions, 
annihilations freeze out when T ~ m x / 20. 

Roughly, the solution to the Boltzmann equation, which tracks the neu¬ 
trino abundance, goes as Y v ~ / ~ ( m{av) an n)~ 1 and hence f \ v h 2 ~ 
(crv) ann 1 , so that parametrically fl„/i 2 ~ 1/m 2 . As a result, the con¬ 
straint on 12 now leads to a lower bound 33,34,35 on the neutrino mass, of 
about m„ >3 — 7 GeV, depending on whether it is a Dirac or Majorana 
neutrino. This bound and the corresponding downward trend ~ 1/m 2 
can again be seen in Figure 4. The result of a more detailed calculation is 
shown in Figure 5 35 for the case of a Dirac neutrino. The two curves show 
the slight sensitivity on the temperature scale associated with the quark- 
hadron transition. The result for a Majorana mass neutrino is qualitatively 
similar. Indeed, any particle with roughly weak scale cross-sections will 
tend to give an interesting value of ilh 2 ~ 1 . 

The deep drop in fl„/i 2 , visible in Figure 4 at around m v = Mz/ 2, is 
due to a very strong annihilation cross section at Z-boson pole. For yet 
higher neutrino masses the Z-annihilation channel cross section drops as 
~ 1/m 2 , leading to a brief period of an increasing trend in fl„/i 2 . However, 
for m„ > mw the cross section regains its parametric form (av) a nn ~ m 2 
due to the opening up of a new annihilation channel to W-boson pairs 36 , 
and the density drops again as fl v h 2 ~ 1/m 2 . The tree level W-channel 
cross section breaks the unitarity at around Cl (few) TeV 37 however, and 
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Figure 5. The relic density of heavy Dirac neutrinos due to annihilations 35 . The curves 
are labeled by the assumed quark-hadron phase transition temperature in MeV. 


the full cross section must be bound by the unitarity limit 38 . This behaves 
again as 1 /m 2 , whereby f \ v h 2 has to start increasing again, until it becomes 
too large again at 200-400 TeV 38 ’ 37 . 

If neutrinos are Dirac particles, and have a nonzero asymmetry the relic 
density could be governed by the asymmetry rather than by the annihilation 
cross section. Indeed, it is easy to see that the neutrino mass density 
corresponding to the asymmetry r\ v = (n v — np)/n 7 is given by 39 p = 
m„77„n 7 , which implies 

n„h 2 ~ 0.004 r/mo (m„/GeV). (16) 

where r] u io = 10 10 7/„. The behaviour of the energy density of neutrinos with 
an asymmetry is shown by the dotted line in the Figure 4. In the figure, we 
have assumed an asymmetry of 7/„ ~ 5 x 10~ n for neutrinos with standard 
weak interaction strength. 

Based on the leptonic and invisible width of the Z boson, experiments at 
LEP have determined that the number of neutrinos is N v = 2.994± 0.0 1 2 40 . 
Conversely, any new physics must fit within these brackets, and thus 
LEP excludes additional neutrinos (with standard weak interactions) with 
masses 777 , „ < 45 GeV. Combined with the limits displayed in Figures 4 and 
5, we see that the mass density of ordinary heavy neutrinos is bound to be 
very small, f 'L u h 2 < 0.001 for masses m v > 45 GeV up to m v ~ 0(100) 
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TeV. Lab constraints for Dirac neutrinos are available 41 , excluding neutri¬ 
nos with masses between 10 GeV and 4.7 TeV. This is significant, since it 
precludes the possibility of neutrino dark matter based on an asymmetry 
between v and v 39 . 

3.3.3. Axions 

Due to space limitations, the discussion of this candidate will be very 
brief. Axions are pseudo-Goldstone bosons which arise in solving the strong 
CP problem 42,43 via a global U(l) Peccei-Quinn symmetry. The invisible 
axion 43 is associated with the flat direction of the spontaneously broken 
PQ symmetry. Because the PQ symmetry is also explicitly broken (the 
CP violating OFF coupling is not PQ invariant) the axion picks up a small 
mass similar to pion picking up a mass when chiral symmetry is broken. 
We can expect that m a ~ fa where f a , the axion decay constant, is 

the vacuum expectation value of the PQ current and can be taken to be 
quite large. If we write the axion field as a = f a 0 , near the minimum, the 
potential produced by QCD instanton effects looks like V ~ ?ti 2 # 2 / 2 . The 
axion equations of motion lead to a relatively stable oscillating solution. 
The energy density stored in the oscillations exceeds the critical density 44 
unless f a < 10 12 GeV. 

Axions may also be emitted stars and supernova 45 . In supernovae, 
axions are produced via nucleon-nucleon bremsstrahlung with a coupling 
CL 4 -/V oc m N /f a - As was noted above the cosmological density limit requires 
fa ^ 10 12 GeV. Axion emission from red giants imply 46 f a > 10 10 GeV 
(though this limit depends on an adjustable axion-electron coupling), the 
supernova limit requires 4 ' / Q > 2 x 10 11 GeV for naive quark model cou¬ 
plings of the axion to nucleons. Thus only a narrow window exists for the 
axion as a viable dark matter candidate. 

4. Supersymmetric Dark Matter 

For the remaining discussion of dark matter, I will restrict my attention to 
supersymmetry and in particular, the minimal supersymmetric standard 
model (MSSM) with R-parity conservation. R-parity is necessary if one 
wants to forbid all new baryon and lepton number violating interactions at 
the weak scale. If R-parity, which distinguishes between “normal” matter 
and the supersymmetric partners and can be defined in terms of baryon, 
lepton and spin as R = (—l) 3B+i+2 ‘ s ', j s unbroken, there is at least one su¬ 
persymmetric particle (the lightest supersymmetric particle or LSP) which 
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must be stable. Thus, the minimal model contains the fewest number of 
new particles and interactions necessary to make a consistent theory. 

There are very strong constraints, however, forbidding the existence of 
stable or long lived particles which are not color and electrically neutral 48 . 
Strong and electromagnetically interacting LSPs would become bound with 
normal matter forming anomalously heavy isotopes. Indeed, there are very 
strong upper limits on the abundances, relative to hydrogen, of nuclear 
isotopes 49 , n/riH < 10 -15 to 10 -29 for 1 GeV < m < 1 TeV. A strongly 
interacting stable relic is expected to have an abundance n/riH ^ 10 -1 ° 
with a higher abundance for charged particles. 

There are relatively few supersymmetric candidates which are not col¬ 
ored and are electrically neutral. The sneutrino 50 is one possibility, but in 
the MSSM, it has been excluded as a dark matter candidate by direct 41 
and indirect 51 searches. In fact, one can set an accelerator based limit on 
the sneutrino mass from neutrino counting, mp > 44.7 GeV 52 . In this 
case, the direct relic searches in underground low-background experiments 
require mp > 20 TeV 41 . Another possibility is the gravitino which is 
probably the most difficult to exclude. I will concentrate on the remaining 
possibility in the MSSM, namely the neutralinos but will return to the case 
of gravitino dark matter as well. 


4.1. Parameters 

The most general version of the MSSM, despite its minimality in particles 
and interactions contains well over a hundred new parameters. The study 
of such a model would be untenable were it not for some (well motivated) 
assumptions. These have to do with the parameters associated with super- 
symmetry breaking. It is often assumed that, at some unification scale, all 
of the gaugino masses receive a common mass, TO 1 / 2 . The gaugino masses 
at the weak scale are determined by running a set of renormalization group 
equations. Similarly, one often assumes that all scalars receive a common 
mass, mo, at the GUT scale. These too are run down to the weak scale. 
The remaining supersymmetry breaking parameters are the trilinear mass 
terms, Aq, which I will also assume are unified at the GUT scale, and the 
bilinear mass term B. There are, in addition, two physical CP violating 
phases which will not be considered here. Finally, there is the Higgs mixing 
mass parameter, /i, and since there are two Higgs doublets in the MSSM, 
there are two vacuum expectation values. One combination of these is re¬ 
lated to the Z mass, and therefore is not a free parameter, while the other 
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combination, the ratio of the two vevs, tan/3, is free. 

The natural boundary conditions at the GUT scale for the MSSM would 
include fi and B in addition to mu 2 , mo, and Aq. In this case, upon running 
the RGEs down to a low energy scale and minimizing the Higgs potential, 
one would predict the values of M z , tan/3 (in addition to all of the sparticle 
masses). Since Mz is known, it is more useful to analyze supersymmetric 
models where Mz is input rather than output. It is also common to treat 
tan /3 as an input parameter. This can be done at the expense of shifting /i 
(up to a sign) and B from inputs to outputs. This model is often referred 
to as the constrained MSSM or CMSSM. Once these parameters are set, 
the entire spectrum of sparticle masses at the weak scale can be calculated. 



Log 10 (Q/GeV) 


Figure 6. RG evolution of the mass parameters in the CMSSM. 

In Fig. 6, an example of the running of the mass parameters in the 
CMSSM is shown. Here, we have chosen mi/ 2 = 250 GeV, mo = 100 GeV, 
tan/3 = 3, A 0 = 0, and p. < 0. Indeed, it is rather amazing that from 
so few input parameters, all of the masses of the supersymmetric particles 
can be determined. The characteristic features that one sees in the figure, 
are for example, that the colored sparticles are typically the heaviest in the 
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spectrum. This is due to the large positive correction to the masses due 
to «3 in the RGE’s. Also, one finds that the B (the partner of the U(l)y 
gauge boson), is typically the lightest sparticle. But most importantly, 
notice that one of the Higgs mass 2 , goes negative triggering electroweak 
symmetry breaking 53 . (The negative sign in the figure refers to the sign of 
the mass 2 , even though it is the mass of the sparticles which are depicted.) 


4.2. Neutralinos 


There are four neutralinos, each of which is a linear combination of the 
R = — 1 neutral fermions 48 : the wino IT 3 , the partner of the 3rd component 
of the SU{2)l gauge boson; the bino, B\ and the two neutral Higgsinos, 
H\ and i? 2 - Assuming gaugino mass universality at the GUT scale, the 
identity and mass of the LSP are determined by the gaugino mass rrii/ 2 , /J, 
and tan/3. In general, neutralinos can be expressed as a linear combination 

X = aB + pW 3 + 7 JT 1 + 6H 2 (17) 


The solution for the coefficients a, (3 ,7 and 6 for neutralinos that make up 
the LSP can be found by diagonalizing the mass matrix 



where Mi (M 2 ) is a soft supersymmetry breaking term giving mass to the 
U(l) (SU(2)) gaugino(s). In a unified theory Mi = M 2 = m i / 2 at the 
unification scale (at the weak scale, Mi ~ |—M 2 ). As one can see, the co¬ 
efficients a, /?, 7 , and S depend only on mp 2 , /i, and tan/3. In the CMSSM, 
the solutions for /r generally lead to a lightest neutralino which is very 
nearly a pure B. 


4.3. The Relic Density 

The relic abundance of LSP’s is determined by solving the Boltzmann 
equation for the LSP number density in an expanding Universe. The 
technique 35 used is similar to that for computing the relic abundance of 
massive neutrinos 33 with the appropriate substitution of the cross section. 
The relic density depends on additional parameters in the MSSM beyond 
mi/ 2 ,/i, and tan/3. These include the sfermion masses, to 7 and the Higgs 
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pseudo-scalar mass, toa, derived from too (and m 1 / 2 )- To determine the 
relic density it is necessary to obtain the general annihilation cross-section 
for neutralinos. In much of the parameter space of interest, the LSP is 
a bino and the annihilation proceeds mainly through sfermion exchange. 
Because of the p-wave suppression associated with Majorana fermions, the 
s-wave part of the annihilation cross-section is suppressed by the outgoing 
fermion masses. This means that it is necessary to expand the cross-section 
to include p-wave corrections which can be expressed as a term proportional 
to the temperature if neutralinos are in equilibrium. Unless the neutralino 
mass happens to lie near near a pole, such as m x ~ mz/2 or mfc/2, in 
which case there are large contributions to the annihilation through di¬ 
rect s-channel resonance exchange, the dominant contribution to the BB 
annihilation cross section comes from crossed t-channel sfermion exchange. 

Annihilations in the early Universe continue until the annihilation rate 
T ~ avn x drops below the expansion rate. The final neutralino relic density 
expressed as a fraction of the critical energy density can be written as 48 


n x h 2 ~ 1.9 x lCT 11 



GeV 

axf + \bx 2 


(19) 


where (T x /T 7 ) 3 accounts for the subsequent reheating of the photon tem¬ 
perature with respect to y, due to the annihilations of particles with mass 
m < Xfm x 54 and Xf = Tf/m x is proportional to the freeze-out temper¬ 
ature. The coefficients a and b are related to the partial wave expansion 
of the cross-section, av = a + bx + .... Eq. (19 ) results in a very good 
approximation to the relic density expect near s-channel annihilation poles, 
thresholds and in regions where the LSP is nearly degenerate with the next 
lightest supersymmetric particle 55 . 


4.4. The CMSSM after WMAP 

For a given value of tan/3, Aq, and sgn(/j,), the resulting regions of accept¬ 
able relic density and which satisfy the phenomenological constraints can 
be displayed on the mi /2 — mo plane. In Fig. 7a, the light shaded region 
corresponds to that portion of the CMSSM plane with tan /3 = 10, Aq = 0, 
and fi > 0 such that the computed relic density yields 0.1 < < 0.3. 

At relatively low values of mi /2 and too, there is a large ‘bulk’ region 
which tapers off as my 2 is increased. At higher values of too, annihilation 
cross sections are too small to maintain an acceptable relic density and 
n x h 2 > 0.3. Although sfermion masses are also enhanced at large my 2 
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(due to RGE running), co-annihilation processes between the LSP and the 
next lightest sparticle (in this case the t 1 ) enhance the annihilation cross 
section and reduce the relic density. This occurs when the LSP and NLSP 
are nearly degenerate in mass. The dark shaded region has < m x and 
is excluded. Neglecting coannihilations, one would find an upper bound of 
~ 450 GeV on my 2 , corresponding to an upper bound of roughly 200 GeV 
on my. The effect of coannihilations is to create an allowed band about 25- 
50 GeV wide in mo for m 1 / 2 ^ 1400 GeV, which tracks above the TOf 1 = m x 
contour 56 . 



Figure 7. The (mi/ 2 ,mo) planes for (a) tan (3 = 10 and p > 0, assuming Aq = 
0, mt = 175 GeV and = 4.25 GeV. The near-vertical (red) dot-dashed lines 

are the contours mh = 114 GeV, and the near-vertical (black) dashed line is the contour 
m x ± = 104 GeV. Also shown by the dot-dashed curve in the lower left is the corner 
excluded by the LEP bound of mg > 99 GeV. The medium (dark green) shaded region is 
excluded by b —► sj, and the light (turquoise) shaded area is the cosmologically preferred 
regions with 0.1 < Ll x h 2 <0.3. In the dark (brick red) shaded region, the LSP is the 
charged fi. The region allowed by the E821 measurement of a^ at the 2-cr level, is 
shaded (pink) and bounded by solid black lines, with dashed lines indicating the 1-a 
ranges. In (b), the relic density is restricted to the range 0.094 < Ll x h 2 < 0.129. 


Also shown in Fig. 7a are the relevant phenomenological constraints. 
These include the limit on the chargino mass: m x ± > 104 GeV 57 , on the 
selectron mass: mg > 99 GeV 58 and on the Higgs mass: rrih >114 GeV 59 . 
The former two constrain mp 2 and mo directly via the sparticle masses, 
and the latter indirectly via the sensitivity of radiative corrections to the 
Higgs mass to the sparticle masses, principally nn FeynHiggs 60 is used 
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for the calculation of rrih. The Higgs limit imposes important constraints 
principally on mp 2 particularly at low tan (3. Another constraint is the 
requirement that the branching ratio for b —> sy is consistent with the ex¬ 
perimental measurements 61 . These measurements agree with the Standard 
Model, and therefore provide bounds on MSSM particles 62,63 , such as the 
chargino and charged Higgs masses, in particular. Typically, the b —» sy 
constraint is more important for g < 0, but it is also relevant for g > 0, 
particularly when tan (3 is large. The constraint imposed by measurements 
of b —> s') also excludes small values of nii/ 2 - Finally, there are regions 
of the (mi /2 ,mo) plane that are favoured by the BNL measurement 64 of 

— 2 at the 2-cr level, corresponding to a deviation from the Standard 
Model calculation 65 using e + e _ data. One should be however aware that 
this constraint is still under active discussion. 

The preferred range of the relic LSP density has been altered signifi¬ 
cantly by the recent improved determination of the allowable range of the 
cold dark matter density obtained by combining WMAP and other cosmo¬ 
logical data: 0.094 < VLqdm < 0.129 at the 2-cr level 9 . In the second panel 
of Fig. 7, we see the effect of imposing the WMAP range on the neutralino 
density 66,67,68 . We see immediately that (i) the cosmological regions are 
generally much narrower, and (ii) the ‘bulk’ regions at small and mo 
have almost disappeared, in particular when the laboratory constraints are 
imposed. Looking more closely at the coannihilation regions, we see that 
(iii) they are significantly truncated as well as becoming much narrower, 
since the reduced upper bound on fl x h 2 moves the tip where m x = mf to 
smaller m- 1/2 so that the upper limit is now TO 1/2 ^ 950 GeV or m x < 400 
GeV. 

Another mechanism for extending the allowed CMSSM region to large 
m x is rapid annihilation via a direct-channel pole when m x ~ \ttia 69,7 °. 
Since the heavy scalar and pseudoscalar Higgs masses decrease as tan (3 
increases, eventually 2 m x ~ iba yielding a ‘funnel’ extending to large mp 2 
and mo at large tan/3, as seen in the high tan/3 strips of Fig. 8. As one 
can see, the impact of the Higgs mass constraint is reduced (relative to the 
case with tan j3 = 10 ) while that of b —■> sj is enhanced. 

Shown in Fig. 9 are the WMAP lines 66 of the (mi/ 2 ,mo) plane allowed 
by the new cosmological constraint 0.094 < fl x h 2 < 0.129 and the labora¬ 
tory constraints listed above, for g > 0 and values of tan /3 from 5 to 55, in 
steps A(tan/3) = 5. We notice immediately that the strips are considerably 
narrower than the spacing between them, though any intermediate point in 
the (m 1( / 2 ,mo) plane would be compatible with some intermediate value of 
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tan B = 50, u > 0 tan B = 50 , li > 0 



100 1000 2000 3000 100 1000 2000 3000 

m 1/2 (GeV) m 1/2 (GeV) 

Figure 8. As in Fig. 7 for tan 0 = 50. 


tan /3. The right (left) ends of the strips correspond to the maximal (mini¬ 
mal) allowed values of mi/ 2 and hence rri x . The lower bounds on mi/ 2 are 
due to the Higgs mass constraint for tan/3 < 23, but are determined by the 
b —> sy constraint for higher values of tan /3. 

Finally, there is one additional region of acceptable relic density known 
as the focus-point region 71 , which is found at very high values of mo- An 
example showing this region is found in Fig. 10, plotted for tan/3 = 10, 
H > 0, and m t = 175 TeV. As mo is increased, the solution for 0 at low 
energies as determined by the electroweak symmetry breaking conditions 
eventually begins to drop. When 0 < m ]/ 2 , the composition of the LSP 
gains a strong Higgsino component and as such the relic density begins to 
drop precipitously. These effects are both shown in Fig. 11 where the value 
of 0 and flh 2 are plotted as a function of mo for fixed mi/ 2 = 300 GeV and 
tan/3 = 10. As mo is increased further, there are no longer any solutions 
for 0. This occurs in the shaded region in the upper left corner of Fig. 10. 

Fig. 11 also exemplifies the degree of fine tuning associated with the 
focus-point region. While the position of the focus-point region in the 
?Bo,mi/ 2 plane is not overly sensitive to supersymmetric parameters, it is 
highly sensitive to the top quark Yukawa coupling which contributes to the 
evolution of 0 72 ’' 3 . As one can see in the figure, a change in m* of 3 GeV 
produces a shift of about 2.5 TeV in mo- Note that the position of the 
focus-point region is also highly sensitive to the value of Aq/tuq. In Fig. 
11, Aq = 0 was chosen. For Ao/mo = 0.5, the focus point shifts from 2.5 
to 4.5 TeV and moves to larger mo as Ao/mo is increased. 
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Figure 9. The strips display the regions of the (mi/ 2 ) m o) plane that are compatible 
with 0.094 < r2 x /i 2 < 0.129 and the laboratory constraints for (i > 0 and tan (3 = 
5,10,15, 20, 25, 30, 35, 40, 45, 50, 55. The parts of the strips compatible with g^ — 2 at the 
2-cr level have darker shading. 


4.5. A Likelihood analysis of the CMSSM 

In displaying acceptable regions of cosmological density in the mo,mi /2 
plane, it has been assumed that the input parameters are known with per¬ 
fect accuracy so that the relic density can be calculated precisely. While all 
of the beyond the standard model parameters are completely unknown and 
therefore carry no formal uncertainties, standard model parameters such as 
the top and bottom Yukawa couplings are known but do carry significant 
uncertainties. 

The optimal way to combine the various constraints (both phenomeno¬ 
logical and cosmological) is via a likelihood analysis. When performing such 
an analysis, in addition to the formal experimental errors, it is also essential 
to take into account theoretical errors, which introduce systematic uncer¬ 
tainties that are frequently non-negligible. Recently, we have preformed an 
extensive likelihood analysis of the CMSSM 74 . 

The interpretation of the combined Higgs likelihood, L exp , in the 
(mi/ 2 ,mo) plane depends on uncertainties in the theoretical calculation 
of rn-ft,. These include the experimental error in mt and (particularly at 
large tan (3) mb, and theoretical uncertainties associated with higher-order 
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m h = 114 GeV 



100 200 300 400 500 600 700 800 900 1000 

m 1/2 (GeV) 


Figure 10. As in Fig. la, where the range in mo is extended to 5 TeV. In the shaded 
region at very high mo, there are no solutions for p which respect the low energy elec- 
troweak symmetry breaking conditions. 

corrections to mh- Our default assumptions are that mt = 175 ±5 GeV for 
the pole mass, and mb = 4.25 ± 0.25 GeV for the running MS mass evalu¬ 
ated at mb itself. The theoretical uncertainty in mh, a t h, is dominated by 
the experimental uncertainties in mt,b, which are treated as uncorrelated 
Gaussian errors: 

Typically, we find that ( dmh./dm t ) ~ 0.5, so that a t h is roughly 2-3 GeV. 

The combined experimental likelihood, C exp , from direct searches at 
LEP 2 and a global electroweak fit is then convolved with a theoretical 
likelihood (taken as a Gaussian) with uncertainty given by ath from (20) 
above. Thus, we define the total Higgs likelihood function, Ch, as 

£h(m h ) = —=— [ dm'h C exp (m' h ) , (21) 

V27 Ta th J 

where J\f is a factor that normalizes the experimental likelihood distribu¬ 
tion. In addition to the Higgs likelihood function, we have included the 
likelihood function based on b —> S'y. While the likelihood function based 
on the measurements of the anomalous magnetic moment of the muon was 
considered in 74 , it will not be discussed here. 
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Figure 11. The value of fi as a function of mo for fixed m-i/2 = 300 GeV and tan/3 = 10 
for two choices of mt as indicated. The scale on the right gives the value of Qh 2 . The 
curves corresponding to this is scale rise sharply at low mo to values much larger than 
1. For mt = 175 GeV and mo ~ 2500 GeV, the value of £lh 2 drops to acceptable values 
when p becomes small. When the mt = 178 GeV, Vlh 2 drops at mo ~ 5000 GeV. 


Finally, in calculating the likelihood of the CDM density, we take into 
account the contribution of the uncertainties in rrit,b • We will see that the 
theoretical uncertainty plays a very significant role in this analysis. The 
likelihood for flh 2 is therefore, 


Cn h 2 


1 -(nh 2th -Clh 2exp ) 2 / 2 a 2 


( 22 ) 


where a 2 = a 2 xp + a 2 h , with a exp taken from the WMAP 9 result and a 2 h 
from (20), replacing m-h by S2 h 2 . 

The total likelihood function is computed by combining all the compo¬ 
nents described above: 


Ctot = Ch x £h S7 x Cn x h 2 (xC ai ,) (23) 

The likelihood function in the CMSSM can be considered a function of two 
variables, £tot{'rni/ 2 i rn o)i where mu 2 an d mo are the unified GUT-scale 
gaugino and scalar masses respectively. 

Using a fully normalized likelihood function £ t ot obtained by combining 
both signs of /r for each value of tan /3, we can determine the regions in the 
( m i/ 2 j m o) planes which correspond to specific CLs as shown in Fig. 12. 
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The darkest (blue), intermediate (red) and lightest (green) shaded regions 
are, respectively, those where the likelihood is above 68%, above 90%, and 
above 95%. 



1000 900 800 700 600 500 400 300 200 100 200 300 400 500 600 700 800 900 1000 

m 1 / 2 (GeV) m 1 / 2 (GeV) 


Figure 12. Contours of the likelihood at the 68%, 90% and 95% levels for tan/3 = 10, 
Ao = 0 and p > 0 (left panel) or p < 0 (right panel), calculated using information of 
nih, b —> S 7 and f 'Icdm h 2 an( t the current uncertainties in mt and m 



m 1 / 2 (GeV) m 1 / 2 (GeV) 


Figure 13. As in Fig. 12 but assuming zero uncertainty in mt. 


The bulk region is more apparent in the right panel of Fig. 12 for fi > 0 
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than it would be if the experimental error in m ( and the theoretical error 
in irih were neglected. Fig. 13 complements the previous figures by showing 
the likelihood functions as they would appear if there were no uncertainty 
in m,t, keeping the other inputs the same. We see that, in this case, both 
the coannihilation and focus-point strips rise above the 68% CL. 


4.6. Beyond the CMSSM 

The results of the CMSSM described in the previous sections are based 
heavily on the assumptions of universality of the supersymmetry breaking 
parameters. One of the simplest generalizations of this model relaxes the 
assumption of universality of the Higgs soft masses and is known as the 
NUHM 75 In this case, the input parameters include y and toa, in addition 
to the standard CMSSM inputs. In order to switch y and toa from outputs 
to inputs, the two soft Higgs masses, mi, m 2 can no longer be set equal 
to mo and instead are calculated from the electroweak symmetry breaking 
conditions. The NUHM parameter space was recently analyzed 75 and a 
sample of the results are shown in Fig. 14. 

In the left panel of Fig. 14, we see a to-j/ 2 ;mo plane with a relative 
low value of y. In this case, an allowed region is found when the LSP 
contains a non-negligible Higgsino component which moderates the relic 
density independent of mo- To the right of this region, the relic density is 
too small. In the right panel, we see an example of the via, y plane. The 
crosses correspond to CMSSM points. In this single pane, we see examples 
of acceptable cosmological regions corresponding to the bulk region, co¬ 
annihilation region and s-channel annihilation through the Higgs pseudo 
scalar. 

Rather than relax the CMSSM, it is in fact possible to further con¬ 
strain the model. While the CMSSM models described above are certainly 
mSUGRA inspired, minimal supergravity models can be argued to be still 
more predictive. In the simplest version of the theory 76 where supersymme¬ 
try is broken in a hidden sector, the universal trilinear soft supersymmetry¬ 
breaking terms are A = (3 — \/3)mo and bilinear soft supersymmetry¬ 
breaking term is B = (2 — v / 3)mo, i.e., a special case of a general relation 
between B and A, Bq = Aq — mo- 

Given a relation between Bq and Aq, we can no longer use the standard 
CMSSM boundary conditions, in which mp 2 , mo, Aq, tan/ 3, and sgn(y) are 
input at the GUT scale with y and B determined by the electroweak sym¬ 
metry breaking condition. Now, one is forced to input Bq and instead tan /3 
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tan B = 35 , u = 400 GeV, in, = 700 GeV tan B = 10 . m,= 300 . m„ = 100 



-2000 -1000 0 

m i/2 (GeV) ji (GeV) 


Figure 14. a) The NUHM (mi/ 2 > m o) plane for tan/3 = 35, (a) pi = 400 GeV and 
mA — 700 GeV b)the NUHM (//, m_a) plane for tan (3 = 10, mo = 100 GeV and m \/2 = 
300 GeV, with Aq = 0. The (red) dot-dashed lines are the contours mh = 114 GeV, 
and the near-vertical (black) dashed lines are the contours m x ± = 103.5 GeV. The dark 
(black) dot-dashed lines indicate the GUT stability constraint. Only the areas inside 
these curves (small fi) are allowed by this constraint. The light (turquoise) shaded 
areas are the cosmologically preferred regions with 0.1 < Ll x h 2 <0.3. The darker (blue) 
portion of this region corresponds to the WMAP densities. The dark (brick red) shaded 
regions is excluded because the stau is the LSP, and the lighter (yellow) shaded regions is 
excluded because the LSP is a sneutrino. The medium (green) shaded region is excluded 
by b —» s'y. The regions allowed by the g — 2 constraint are shaded (pink) and bounded 
by solid black lines. The solid (blue) curves correspond to m x = 7714 / 2 . 

is calculated from the minimization of the Higgs potential 77 . In Fig. 15, the 
contours of tan/3 (solid blue lines) in the (toi/ 2 ,too) planes for two values 
of A = A 0 /m o, B = Bo/mo = A —1 and the sign of /i are displayed 7 '. 

In panel (a) of Fig. 15, we see that the Higgs constraint combined with 
the relic density requires tan (3 > 11, whilst the relic density also enforces 
tan/3 < 20. For a given point in the toi /2 — mo plane, the calculated value 
of tan /3 increases as A increases. This is seen in panel (b) of Fig. 15, when 
A = 2.0, close to its maximal value for /i > 0, the tan/3 contours turn over 
towards smaller mi/ 2 , and only relatively large values 25 < tan/3 < 35 
are allowed by the b —> S'y and flcDMh 2 constraints, respectively. For any 
given value of A , there is only a relatively narrow range allowed for tan/3. 

4.7. Detectability 

Direct detection techniques rely on an ample neutralino-nucleon scatter¬ 
ing cross-section. In Fig. 16, we display the allowed ranges of the spin- 
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Figure 15. Examples of (mi/ 2 , mo) planes with contours o/tan/3 superposed, for p > 0 
and (a) the simplest Polonyi model with A = 3 — v3, B = A — 1 and (b) A = 2.0, B = 
A — 1. /n each panel, we show the regions excluded by the LEP lower limits on MSSM 
particles, those ruled out by b —> s'y decay (medium green shading), and those excluded 
because the LSP would be charged (dark red shading). The region favoured by the WMAP 
range has light turquoise shading. The region suggested by — 2 is medium (pink) 
shaded. 

independent cross sections in the NUHM when we sample randomly tan/? 
as well as the other NUHM parameters 78 . The raggedness of the boundaries 
of the shaded regions reflects the finite sample size. The dark shaded regions 
includes all sample points after the constraints discussed above (including 
the relic density constraint) have been applied. In a random sample, one 
often hits points which are are perfectly acceptable at low energy scales 
but when the parameters are run to high energies approaching the GUT 
scale, one or several of the sparticles mass squared runs negative. This has 
been referred to as the GUT constraint here. The medium shaded region 
embodies those points after the GUT constraint has been applied. After 
incorporating all the cuts, including that motivated by — 2, we find that 
the light shaded region where the scalar cross section has the range 10 -6 pb 
> &si 10 -10 pb, with somewhat larger (smaller) values being possible in 
exceptional cases. 

The results from this analysis 78 for the scattering cross section in the 
NUHM (which by definition includes all CMSSM results) are compared 
with the previous CDMS 79 and Edelweiss 80 bounds as well as the recent 
CDMSII results 81 in Fig. 16. While previous experimental sensitivities were 
not strong enough to probe predictions of the NUHM, the current CDMSII 
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Figure 16. Ranges of the spin-independent cross section in the NUHM. The ranges 
allowed by the cuts on D x h 2 , m ^ and b —*■ S 7 have dark shading, those still allowed 
by the GUT stability cut have medium shading, and those still allowed after applying 
all the cuts including g^ — 2 have light shading. The pale shaded region corresponds to 
the extra area of points with low relic densities, whose cross sections have been rescaled 
appropriately. Also shown are the limits from the CDMS 79 and Edelweiss 80 experiments 
as well as the recent CDMSII result 81 on the neutralino-proton elastic scattering cross 
section as a function of the neutralino mass. The CDMSII limit is stronger than the 
Edelweiss limit which is stronger than the previous CDMS limit at higher m x . The 
result reported by DAM A 82 is found in the upper left. 


bound has begun to exclude realistic models and it is expected that these 
bounds improve by a factor of about 20. See ref. 83 for updated direct 
detection calculations in the MSSM. 

5. Big Bang Nucleosynthesis 

The standard model 84 of big bang nucleosynthesis (BBN) is based on the 
relatively simple idea of including an extended nuclear network into a ho¬ 
mogeneous and isotropic cosmology. Apart from the input nuclear cross 
sections, the theory contains only a single parameter, namely the baryon- 
to-photon ratio, rj. Other factors, such as the uncertainties in reaction 
rates, and the neutron mean-life can be treated by standard statistical and 
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Monte Carlo techniques 85,86,87,88,89 . The theory then allows one to make 
predictions (with well-defined uncertainties) of the abundances of the light 
elements, D, 3 He, 4 He, and 7 Li. 

5.1. Theory 

Conditions for the synthesis of the light elements were attained in the early 
Universe at temperatures T > 1 MeV. In the early Universe, the energy 
density was dominated by radiation with 

' , = 55 (2 + U3 jv - )t4 < 24 > 

from the contributions of photons, electrons and positrons, and N„ neutrino 
flavors (at higher temperatures, other particle degrees of freedom should be 
included as well). At these temperatures, weak interaction rates were in 
equilibrium. In particular, the processes 

n + e + <-> p + v e 
n + v e <-> p + e~ 

n <-> p + e~ + u e (25) 

fix the ratio of number densities of neutrons to protons. At T 1 MeV, 

(n/p) ~ 1. 

The weak interactions do not remain in equilibrium at lower temper¬ 
atures. Freeze-out occurs when the weak interaction rate, T^,*, ~ G 2 F T b 
falls below the expansion rate which is given by the Hubble parameter, 
H ~ VGFp ~ T 2 /Mp, where M P = 1/VGn - 1.2 x 10 19 GeV. The 
/3-interactions in eq. (25) freeze-out at about 0.8 MeV. As the tempera¬ 
ture falls and approaches the point where the weak interaction rates are 
no longer fast enough to maintain equilibrium, the neutron to proton ratio 
is given approximately by the Boltzmann factor, (n/p) ~ e _Am / T ~ 1/6, 
where A?n is the neutron-proton mass difference. After freeze-out, free 
neutron decays drop the ratio slightly to about 1/7 before nucleosynthesis 
begins. 

The nucleosynthesis chain begins with the formation of deuterium by 
the process, p + n — ■> D + 7. However, because of the large number of 
photons relative to nucleons, = n 7 /ns ~ 10 10 , deuterium production 
is delayed past the point where the temperature has fallen below the deu¬ 
terium binding energy, E p = 2.2 MeV (the average photon energy in a 
blackbody is £/, ~ 2.7 T). This is because there are many photons in the 
exponential tail of the photon energy distribution with energies E > Ep 
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Figure 17. The light element abundances from big bang nucleosynthesis as a function 
of 1 ) 10 . 


despite the fact that the temperature or E 1 is less than Eg- The degree 
to which deuterium production is delayed can be found by comparing the 
qualitative expressions for the deuterium production and destruction rates, 


T p ss ngcrv (26) 

ss n 1 ave~ EB ^ T 


When the quantity ?y _1 exp(— Eg/T) ~ 1, the rate for deuterium destruc¬ 
tion (D + 7 —> p + n) finally falls below the deuterium production rate and 
the nuclear chain begins at a temperature T ~ 0.1 MeV. 

The dominant product of big bang nucleosynthesis is 4 He and its abun¬ 
dance is very sensitive to the (n/p) ratio 


^ = 2 (n/p) 

p [i + Np)\ 


0.25 


(27) 


i.e., an abundance of close to 25% by mass. Lesser amounts of the other 
light elements are produced: D and 3 He at the level of about 10 -5 by 
number, and 'Li at the level of 10~ 10 by number. 
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Recently the input nuclear data 

have been carefully reassessed 87,88,89,90,91 , leading to improved precision 
in the abundance predictions. The NACRE collaboration presented an up¬ 
dated nuclear compilation 90 . For example, notable improvements include a 
reduction in the uncertainty in the rate for 3 He(?r,p)T from 10% 86 to 3.5% 
and for T(a, 7 ) 7 Li from ~ 23 — 30% 86 to ~ 4%. Since then, new data and 
techniques have become available, motivating new compilations. Within 
the last year, several new BBN compilations have been presented 92,93,94 . 

The resulting elemental abundances predicted by standard BBN are 
shown in Fig. 17 as a function of y 88 . The left plot shows the abundance 
of 4 He by mass, Y, and the abundances of the other three isotopes by 
number. The curves indicate the central predictions from BBN, while the 
bands correspond to the uncertainty in the predicted abundances. This 
theoretical uncertainty is shown explicitly in the right panel as a function 
of T]. 

In the standard model with N„ = 3, the only free parameter is the 
density of baryons which sets the rates of the strong reactions. Thus, any 
abundance measurement determines ??, while additional measurements over¬ 
constrain the theory and thereby provide a consistency check. BBN has 
thus historically been the premier means of determining the cosmic baryon 
density. With the increased precision of microwave background anisotropy 
measurements, it is now possible to use the the CMB to independently 
determine the baryon density. The WMAP value for fish 2 = 0.0224 trans¬ 
lates into 

7710 = 6.14 ± 0.25 (28) 

With 77 fixed by the CMB, precision comparisons to the observations can 
now be attempted 95 . 

5.2. Light Element Observations and Comparison with 
Theory 

BBN theory predicts the universal abundances of D, 3 He, 4 He, and 7 Li, 
which are essentially determined by t ~ 180 s. Abundances are however 
observed at much later epochs, after stellar nucleosynthesis has commenced. 
The ejected remains of this stellar processing can alter the light element 
abundances from their primordial values, and produce heavy elements such 
as C, N, O, and Fe (“metals”). Thus one seeks astrophysical sites with low 
metal abundances, in order to measure light element abundances which 
are closer to primordial. For all of the light elements, systematic errors 
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are an important and often dominant limitation to the precision of derived 
primordial abundances. 

5.2.1. D/H 

In recent years, high-resolution spectra have revealed the presence of D 
in high-redshift, low-metallicity quasar absorption systems (QAS), via its 
isotope-shifted Lyman-a absorption. These are the first measurements of 
light element abundances at cosmological distances. It is believed that there 
are no astrophysical sources of deuterium 96 , so any measurement of D/H 
provides a lower limit to primordial D/H and thus an upper limit on 77 ; 
for example, the local interstellar value of D/H=(1.5 ± 0.1) x 10 -5 citelin 
requires that 7710 < 9. In fact, local interstellar D may have been depleted 
by a factor of 2 or more due to stellar processing; however, for the high- 
redshift systems, conventional models of galactic nucleosynthesis (chemical 
evolution) do not predict significant D/H depletion 98 . 

The five most precise observations of deuterium 99,100,101,102 in QAS give 
D/H = (2.78 ± 0.29) x 10 -5 , where the error is statistical only. These are 
shown in Fig. 18 along with some other recent measurements 103,104,105 . 
Inspection of the data shown in the figure clearly indicates the need for 
concern over systematic errors. We thus conservatively bracket the observed 
values with a range D/H = 2 — 5 x 10 -5 which corresponds to a range in 
7710 of 4 - 8 which easily brackets the CMB determined value. 

Using the WMAP value for the baryon density (28) the primordial D/H 
abundance is predicted to be 88,92 : 

(D/H) p = 2.55tg;|) x 10 " 5 ( 29 ) 

As one can see, this value is in very good agreement with the observational 
value. 

5.3. 4 He 

We observe 4 He in clouds of ionized hydrogen (HII regions), the most metal- 
poor of which are in dwarf galaxies. There is now a large body of data on 
4 He and CNO in these systems 106 . Of the modern 4 He determinations, the 
work of Pagel et al. 10 ' established the analysis techniques that were soon 
to follow 108 . Their value of Y p = 0.228 ± 0.005 was significantly lower than 
that of a sample of 45 low metallicity HII regions, observed and analyzed 
in a uniform manner 106 , with a derived value of Y p = 0.244 ± 0.002. An 
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Figure 18. D/H abundances shown as a function of [Si/H], 


analysis based on the combined available data as well as unpublished data 
yielded an intermediate value of 0.238 ± 0.002 with an estimated systematic 
uncertainty of 0.005 109 . An extended data set including 89 HII regions 
obtained Y p = 0.2429 ± 0.0009 110 . However, the recommended value is 
based on the much smaller subset of 7 HII regions, finding Y p = 0.2421 ± 
0 . 0021 . 

4 He abundance determinations depend on a number of physical parame¬ 
ters associated with the HII region in addition to the overall intensity of the 
He emission line. These include, the temperature, electron density, optical 
depth and degree of underlying absorption. A self-consistent analysis may 
use multiple 4 He emission lines to determine the He abundance, the electron 
density and the optical depth. In 106 , five He lines were used, underlying 
He absorption was assumed to be negligible and used temperatures based 
on OIII data. 

The question of systematic uncertainties was addressed in some detail 
in 111 . It was shown that there exist severe degeneracies inherent in the 
self-consistent method, particularly when the effects of underlying absorp- 
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tion are taken into account. The results of a Monte-Carlo reanalysis 112 of 
NCG 346 113 is shown in Fig. 19. In the left panel, solutions for the 4 He 
abundance and electron density are shown (symbols are described in the 
caption). In the right panel, a similar plot with the 4 He abundance and 
the equivalent width for underlying absorption is shown. As one can see, 
solutions with no absorption and high density are often indistinguishable 
(i.e., in a statistical sense they are equally well represented by the data) 
from solutions with underlying absorption and a lower density. In the lat¬ 
ter case, the He abundance is systematically higher. These degeneracies 
are markedly apparent when the data is analyzed using Monte-Carlo meth¬ 
ods which generate statistically viable representations of the observations 
as shown in Fig. 19. When this is done, not only are the He abundances 
found to be higher, but the uncertainties are also found to be significantly 
larger than in a direct self-consistent approach. 
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Figure 19. Results of modeling of 6 He I line observations of NGC 346 113 . The solid 
lines show the original derived values and the dashed lines show the 1 cr errors on those 
values. The solid circles (with error bars) show the results of the % 2 minimization 
solution (with calculated errors) 112 . The small points show the results of Monte Carlo 
realizations of the original input spectrum. The solid squares (with error bars) show 
the means and dispersions of the output values for the x 2 minimization solutions of the 
Monte Carlo realizations. 


Recently a careful study of the systematic uncertainties in 4 He, particu¬ 
larly the role of underlying absorption has been performed using a subset of 
the highest quality from the data of Izotov and Thuan 106 . All of the phys¬ 
ical parameters listed above including the 4 He abundance were determined 
self-consistently with Monte Carlo methods 111 . Note that the 4 He abun¬ 
dances are systematically higher, and the uncertainties are several times 
larger than quoted in 106 . In fact this study has shown that the determined 























February 2, 2008 16:44 Proceedings Trim Size: 9in x 6in 


otasi4up 


34 

value of Y p is highly sensitive to the method of analysis used. The result is 
shown in Fig. 20 together with a comparison of the previous result. The 
extrapolated 4 He abundance was determined to be Y p = 0.2495 ± 0.0092. 
The value of 77 corresponding to this abundance is 7710 = 6.91 4 1 q and clearly 
overlaps with t]cmb- Conservatively, it would be difficult at this time to 
exclude any value of Y p inside the range 0.232 - 0.258. 

At the WMAP value for 77 , the 4 He abundance is predicted to be 88,92 : 

Y p = 0.2485 ± 0.0005 (30) 

This value is considerably higher than any prior determination of the pri¬ 
mordial 4 He abundance, it is in excellent agreement with the most recent 
analysis of the 4 He abundance 112 . Note also that the large uncertainty as¬ 
cribed to this value indicates that the while 4 He is certainly consistent with 
the WMAP determination of the baryon density, it does not provide for a 
highly discriminatory test of the theory at this time. 



Figure 20. A comparison of the results for the best targets 106 and a re-analysis of the 
spectra for those targets 112 . 


5.4. 7 Li/H 

The systems best suited for Li observations are metal-poor halo stars in our 
Galaxy. Observations have long shown 114 that Li does not vary significantly 
in Pop II stars with metallicities < 1/30 of solar — the “Spite plateau”. 
Recent precision data suggest a small but significant correlation between 
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Li and Fe 115 which can be understood as the result of Li production from 
Galactic cosmic rays 116 . Extrapolating to zero metallicity one arrives at a 
primordial value 117 Li/H| p = (1.23±o; 8 g) x 10~ 10 . 



Figure 21. Contributions to the total predicted lithium abundance from the adopted 
GCE model of 118 , compared with low metallicity stars and a sample of high metallicity 
stars. The solid curve is the sum of all components. 


Figure 21 shows the different Li components for a model with 
( 7 Li/H) p = 1.23 x 10 _1 °. The linear slope produced by the model is in¬ 
dependent of the input primordial value. The model of ref. 118 includes 
in addition to primordial 7 Li, lithium produced in galactic cosmic ray nu¬ 
cleosynthesis (primarily a + a fusion), and ' Li produced by the ^-process 
during type II supernovae. As one can see, these processes are not sufficient 
to reproduce the population I abundance of 'Li, and additional production 
sources are needed. 

Recent data 119 with temperatures based on Ha lines (considered to give 
systematically high temperatures) yields 7 Li/H = (2.19 ± 0.28) x 10 _1 °. 
These results are based on a globular cluster sample (NGC 6397). This 
result is consistent with previous Li measurements of the same cluster which 
gave 7 Li/H = (1.91T0.44) x 10~ 10 120 and 7 Li/H = (1.69T0.27) x 10" 10 121 . 
A related study (also of globular cluster stars) gives ' Li/H = (2.29±0.94) x 

q —10 122 

The 7 Li abundance based on the WMAP baryon density is predicted to 
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be 88 ’ 92 : 


7 Li/H = 4.26 toil x 10_1 ° ( 31 ) 

This value is in clear contradiction with most estimates of the primordial 
Li abundance. It is a factor of ~ 3 higher than the value observed in most 
halo stars, and just about 0.2 dex over the globular cluster value. 


5.5. Concordance 

In Fig. 22, we show the direct comparison between the BBN predicted 
abundances given in eqs. (29), (30), and (31), using the WMAP value of 
rjio = 6.25 ± 0.25 with the observations 123 . As one can see, there is very 
good agreement between theory and observation for both D/H and 4 He. Of 
course, in the case of 4 He, concordance is almost guaranteed by the large 
errors associated to the observed abundance. In contrast, as was just noted 
above, there is a marked discrepancy in the case of 7 Li. 



Figure 22. Primordial light element abundances as predicted by BBN and WMAP 
(dark shaded regions) 123 . Different observational assessments of primordial abundances 
are plotted as follows: (a) the light shaded region shows D/H = (2.78 ± 0.29) X 10 -5 ; 
(b) the light shaded region shows Y p = 0.249 ± 0.009; (c) the light shaded region shows 
7 Li/H = 1.23j)Q' 3 g X 10 _1 °, while the dashed curve shows 7 Li/H = (2.19T0.28) X 10 -10 . 


The quoted value for the 7 Li abundance assumes that the Li abundance 
in the stellar sample reflects the initial abundance at the birth of the star. 
However, an important source of systematic uncertainty comes from the 
possible depletion of Li over the > 10 Gyr age of the Pop II stars. The 
atmospheric Li abundance will suffer depletion if the outer layers of the 
stars have been transported deep enough into the interior, and/or mixed 
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with material from the hot interior; this may occur due to convection, 
rotational mixing, or diffusion. Standard stellar evolution models predict 
Li depletion factors which are very small (<0.05 dex) in very metal-poor 
turnoff stars 124 . However, there is no reason to believe that such simple 
models incorporate all effects which lead to depletion such as rotationally- 
induced mixing and/or diffusion. Current estimates for possible depletion 
factors are in the range ~ 0.2-0.4 dex 125 . As noted above, this data 
sample 115 shows a negligible intrinsic spread in Li leading to the conclusion 
that depletion in these stars is as low as 0.1 dex. 

Another important source for potential systematic uncertainty stems 
from the fact that the Li abundance is not directly observed but rather, 
inferred from an absorption line strength and a model stellar atmosphere. 
Its determination depends on a set of physical parameters and a model- 
dependent analysis of a stellar spectrum. Among these parameters, are the 
metallicity characterized by the iron abundance (though this is a small ef¬ 
fect), the surface gravity which for hot stars can lead to an underestimate of 
up to 0.09 dex if log g is overestimated by 0.5, though this effect is negligible 
in cooler stars. Typical uncertainties in log g are ±0.1 — 0.3. The most im¬ 
portant source for error is the surface temperature. Effective-temperature 
calibrations for stellar atmospheres can differ by up to 150-200 K, with 
higher temperatures resulting in estimated Li abundances which are higher 
by ~ 0.08 dex per 100 K. Thus accounting for a difference of 0.5 dex be¬ 
tween BBN and the observations, would require a serious offset of the stellar 
parameters. While there has been a recent analysis 126 which does support 
higher temperatures, the consequences of the higher temperatures on the 
inferred abundances of related elements such as Be, B, and O observed in 
the same stars is somewhat negative 127 . 

Finally a potential source for systematic uncertainty lies in the BBN 
calculation of the ' Li abundance. As one can see from Fig. 17, the pre¬ 
dictions for ' Li carry the largest uncertainty of the 4 light elements which 
stem from uncertainties in the nuclear rates. The effect of changing the 
yields of certain BBN reactions was recently considered by Coc et al. 91 . 
In particular, they concentrated on the set of cross sections which affect 
' Li and are poorly determined both experimentally and theoretically. In 
many cases however, the required change in cross section far exceeded any 
reasonable uncertainty. Nevertheless, it may be possible that certain cross 
sections have been poorly determined. In 91 , it was found for example, 
that an increase of either the ' Li(d, n)2 4 He or 7 Be(d,p)2 4 He reactions by a 
factor of 100 would reduce the 7 Li abundance by a factor of about 3. 
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The possibility of systematic errors in the 3 He(a, 7 ) 'Be reaction, which 
is the only important 'Li production channel in BBN, was considered in 
detail in 128 . The absolute value of the cross section for this key reaction 
is known relatively poorly both experimentally and theoretically. How¬ 
ever, the agreement between the standard solar model and solar neutrino 
data thus provides additional constraints on variations in this cross section. 
Using the standard solar model of Bahcall 129 , and recent solar neutrino 
data 130 , one can exclude systematic variations of the magnitude needed to 
resolve the BBN 'Li problem at the > 95% CL 128 . Thus the “nuclear fix” 
to the 7 Li BBN problem is unlikely. 

Finally, we turn to 3 He. Here, the only observations available are 
in the solar system and (high-metallicity) HII regions in our Galaxy 131 . 
This makes inference of the primordial abundance difficult, a problem com¬ 
pounded by the fact that stellar nucleosynthesis models for 3 He are in con¬ 
flict with observations 132 . Consequently, it is not appropriate to use 3 He 
as a cosmological probe 133 ; instead, one might hope to turn the problem 
around and constrain stellar astrophysics using the predicted primordial 
3 He abundance 134 . For completeness, we note that the 3 He abundance is 
predicted to be: 

3 He/H = 9.28±g;|t x 10 ” 6 (32) 

at the WMAP value of 77 . 


6. Constraints on Decaying Particles and Gravitino Dark 

Matter from BBN 

As an example of constraints on particle properties from BBN, I will con¬ 
centrate here on life-time and abundance limits on decaying particles as it 
ties in well with the previous discussion on supersymmetric dark matter. 
There are of course many other constraints on particle properties which can 
be derived from BBN, most notably the limit on the number of relativistic 
degrees of freedom. For a recent update on these limits, see 123 . 

Because there is good overall agreement between the theoretical pre¬ 
dictions of the light element abundances and their observational determi¬ 
nation, any departure from the standard model (or either particle physics, 
cosmology, or BBN) generally leads to serious inconsistencies among the 
element abundances. 

Gravitinos have long been known to be potentially problematic in 
cosmology 135 . If gravitinos are present with equilibrium number densities, 
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we can write their energy density as 


P3/2 — m 3/2 n 3/2 


m 3 / 2 



(33) 


where today one expects that the gravitino temperature T 3 / 2 is reduced 
relative to the photon temperature due to the annihilations of particles 
dating back to the Planck time 54 . Typically one can expect the gravitino 
abundance Y 3 / 2 = n 3 / 2 /n 1 ~ (T 3 / 2 /T 7 ) 3 ~ 10” 2 . Then for fl 3 / 2 h 2 < 1, we 
obtain the limit that m 3 / 2 < 1 keV. 

Of course, the above mass limit assumes a stable gravitino, the problem 
persists however, even if the gravitino decays, since its gravitational decay 
rate is very slow. Gravitinos decay when their decay rate, T 3 / 2 ~ m|/ 2 /Mj>, 
becomes comparable to the expansion rate of the Universe (which becomes 
dominated by the mass density of gravitinos), H ~ Thus 

decays occur at T d ~ m^/Mp 3 . After the decay, the Universe is “re¬ 
heated” to a temperature 

T r ~ p(T d y^ ~ rriyz/Mp 2 (34) 

The Universe must reheat sufficiently so that big bang nucleosynthesis oc¬ 
curs in a standard radiation dominated Universe. For Tr > 1 MeV, we 
must require in 3 / 2 > 20 TeV. This large value threatens the solution of the 
hierarchy problem. 

Inflation could alleviate the gravitino problem by diluting the gravitino 
abundance to safe levels 136 . If gravitinos satisfy the noninflationary bounds, 
then their reproduction after inflation is never a problem. For gravitinos 
with mass of order 100 GeV, dilution without over-regeneration will also 
solve the problem, but there are several factors one must contend with in 
order to be cosmologically safe. Gravitino decay products can also upset 
the successful predictions of Big Bang nucleosynthesis, and decays into 
LSPs (if R-parity is conserved) can also yield too large a mass density in 
the now-decoupled LSPs 48 . For unstable gravitinos, the most restrictive 
bound on their number density comes form the photo-destruction of the 
light elements produced during nucleosynthesis 137,138,139 . 

Here, we will consider electromagnetic decays, meaning that the decays 
inject electromagnetic radiation into the early universe. If the decaying par¬ 
ticle is abundant enough or massive enough, the injection of electromagnetic 
radiation can photo-erode the light elements created during primordial nu¬ 
cleosynthesis. The theories we have in mind are generally supersymmetric, 
in which the gravitino and neutralino are the next-to-lightest and lightest 
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supersymmetric particles, respectively (or vice versa), but the constraints 
hold for any decay producing electromagnetic radiation. We thus constrain 
the abundance of such a particle given its mean lifetime tx ■ The abundance 
is constrained through the parameter (/x = mx n x/n-y 

The BBN limits in the (Qx ■, tx) plane is shown in Fig. 23 140 . The 
constraint placed by the 4 He abundance comes from its lower limit, as this 
scenario destroys 4 He. Shown are the limits assuming Y m i n = 0.232 and 
0.22 7 138 » 1 23,i4o^ Tbe area a p, ove these curves are excluded. The deuterium 
lines correspond to the contours (1.3 or 2.2) x 10” 5 < D/H < 5.3xl0“ 5 . 
The first of the lower bounds is the higher line to the left of the cleft, and 
represents the very conservative lower limit on D/H assumed in 138 . The 
range 1.3 - 5.3 xlO -5 effectively brackets all recent observations of D/H 
in quasar absorption systems as discussed above. The second of the lower 
bounds is the lower line on the left side and represents the 2-tr lower limit 
in the best set of D/H observations. The upper bound is the line to the 
right of the cleft. A priori , there is also a narrow strip at larger (,x and 
tx where the D/H ratio also falls within the acceptable range but this is 
excluded by the observed 4 He abundance. 



Figure 23. The constraints imposed by the astrophysical observations of 4 He (red lines), 
D/H (green lines), 6 Li (yellow line), 6 Li/ 7 Li (blue lines), ‘Li (blue band) and 3 He 
(black lines). 
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The constraint imposed by the 6 Li abundance is shown 138 as a solid 
yellow line in Fig. 23. Also shown, as solid blue lines, are two contours 
representing upper limits on the 6 Li/ 7 Li ratio: 6 Li/ 7 Li < 0.07 or 0.15. 
The lower number was used in 138 and represented the upper limit available 
at the time, which was essentially based on multiple observations of a single 
star. The most recent data 141 includes observations of several stars. The Li 
isotope ratio for most metal-poor stars in the sample is as high as 0.15, and 
we display that upper limit here 140 . The main effect of this constraint is 
to disallow a region in the near-vertical cleft between the upper and lower 
limits on D/H, as seen in Fig. 23. 

The blue shaded band in Fig. 23 corresponds to a 7 Li abundance of 
0.9 x 10 -10 < 7 Li/H < (2 or 3) x 10 -10 , with the 7 Li abundance 

decreasing as Qx increases and the intensity of the shading changing at the 
intermediate value. In 138 , only the lower bound was used due the existing 
discrepancy between the primordial and observationally determined values. 
It is apparent that 7 Li abundances in the lower part of the range are possible 
only high in the Deuterium cleft, and even then only if one uses the recent 
and more relaxed limit on the 6 Li/ 7 Li ratio. Values of the 7 Li abundance 
in the upper part of the range are possible, however, even if one uses the 
more stringent constraint on 6 Li/ 7 Li. In this case, the allowed region of 
parameter space would also extend to lower tx , if one could tolerate values 
of D/H between 1.3 and 2.2 x 10 -5 . 

Finally, we show the impact of the 3 He constraint 139,140 . Since Deu¬ 
terium is more fragile than 3 He, whose abundance is thought to have re¬ 
mained roughly constant since primordial nucleosynthesis when comparing 
the BBN value to it proto-solar abundance, one would expect, in principle, 
the 3 He/D ratio to have been increased by stellar processing. Since D is 
totally destroyed in stars, the ratio of 3 He/D can only increase in time or 
remain constant if 3 He is also completely destroyed in stars. The present 
or proto-solar value of 3 He/D can therefore be used to set an upper limit 
on the primordial value. Fig. 23 displays the upper limits 3 He/Z? < 1 or 2 
as solid black lines. Above these contours, the value of 3 He/D increases 
very rapidly, and points high in the Deuterium cleft of Fig. 23 have ab¬ 
surdly high values of 3 He/D, exceeding the limit by an order of magnitude 
or more. 

The previous upper limit on rjx 138 corresponded to the constraint 
mxnx/n-y < 5.0 x 10 -12 GeV for tx = 10 8 s. The weaker (stronger) 
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version of the 3 He constraint adopted corresponds 140 to 

m x — < 2.0(0.8) x 10~ 12 GeV (35) 

n 7 

for rx = 10 8 s. 

Returning to the case of a decaying gravitino, recall that thermal reac¬ 
tions are estimated to produce an abundance of gravitinos given by 142,138 : 

^ = (O.T-2.7)xl0-"x( T5 2i^). (36) 

Assuming that ?n 3 / 2 = 100 GeV and tx = 10® s, and imposing the con¬ 
straints (35), we now find 

T r < (0.8 - 2.8) x 10 7 GeV, ((0.3 - 1.1) x 10 7 GeV) (37) 

for the weaker (stronger) version of the 3 He constraint. 

Finally, we consider the possibility that gravitinos are stable and the 
LSP 143,144 . In this case, in the CMSSM, the next lightest supersymmetric 
particle (NSP) is either the neutralino or the stau. In Fig. 24, we fix the 
ratio of supersymmetric Higgs vacuum expectation values tan (3 = 10 (left 
panel), and tan/3 = 57 (right panel), and assume m 3 / 2 = 100 GeV. In each 
panel of Fig. 24, we display accelerator, astrophysical and cosmological 
constraints in the corresponding (mi/ 2 ,mo) planes as discussed above for 
the CMSSM, concentrating on the regions to the right of the near-vertical 
black lines, where the gravitino is the LSP. The NSP is the r lepton below 
the (red) dotted line. 

Below and to the right of the upper (purple) dashed lines, the density 
of relic gravitinos produced in the decays of other supersymmetric parti¬ 
cles is always below the WMAP upper limit: fl 3 / 2 /i 2 < 0.129. The code 
used in 138 , when combined with the observational constraints used in 138 , 
yielded the astrophysical constraint represented by the dashed grey-green 
lines in both panels of Fig. 24 and did not include the constraint due to 
3 He/D. These constraints on the CMSSM parameter plane were computed 
in 144 . For each point in the (mi/ 2 ,mo), the relic density of either y or f 
is computed and Qx is determined using -h 2 = 3.9 x 10 7 GeV Qx- When 
X = f, (x is reduced by a factor of 0.3, as only 30% of stau decays result 
in electromagnetic showers which affect the element abundances at these 
lifetimes. In addition, at each point, the lifetime of the NSP is computed. 
Then for each t\, the limit on £x is found from the results shown in Fig. 
23. The region to the right of this curve where r = C,x /Cx mlt < 1 is allowed. 
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Figure 24. The (mi/ 2 ,mo) planes for p > 0, 7713/2 = 100 GeV and (a) tan (3 = 10 (b) 
tan (3 = 57. We restrict our attention to the regions between the solid black lines, where 
the gravitino is the LSP and the NSP lifetime exceeds 10 4 s. In each panel, the near¬ 
vertical dashed black (dash-dotted red) line is the constraint m x ± > 104 GeV (m^ > 
114 GeV), the upper (purple) dashed line is the constraint Vl^^h 2 < 0.129, and the light 
green shaded region is that where the NSP would have had 0.094 < flh 2 < 0.129 if it 
had not decayed. The solid red (dashed grey-green) line is the region now (previously) 
allowed by the light-element abundances: r < 1 as described in the text. The red (blue) 
shaded region is that where the 7 Li abundance could have been improved by NSP decays, 
but which is now excluded by the 3 He (D) constraint. 


The astrophysical constraints obtained with the newer abundance 
limits 140 yields the solid red lines in Fig. 24. The examples where tx 
and (x for the NSP decays fall within the ranges shown by the blue band 
of Fig. 23, and hence are suitable for modifying the 'Li abundance 145,146 , 
are shown as red and blue shaded regions in each panel of Fig. 24. If we had 
been able to allow a Deuterium abundance as low as D/H ~ (1 — 2) x 10 -5 , 
the blue shaded region would have been able to resolve the Li discrepancy 
in the context of the CMSSM with gravitino dark matter. The blue region 
that we now regard as excluded by the lower limit on D/H, which is stronger 
than that used in 138 , extends to large m 1 / 2 . The red shaded region, which 
is consistent even with this limit on D/H, but yields very large 3 He/D. 
Fig. 24 show as solid red lines the additional restrictions these constraints 
impose on the (mirrto) planes 140 . 


7. The variation of fundamental constants 

There has been considerable interest of late in the possible variation of the 
fundamental constants. The construction of theories with variable “con- 
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stants” is straightforward. Consider for example a gravitational Lagrangian 
which contains the term 


C ~ <j>R, 


(38) 


where </> is some scalar field and R is the Einstein curvature scalar. The 
gravitational constant is determined if the dynamics of the theory fix the 
expectation value of the scalar field so that 


G n = 


167r(</>) 


(39) 


Similarly a coupling in the Lagrangian of a scalar to the Maxwell term F 2 , 
fixes the fine-structure constant 


C~<t>F 2 , 


1 

167r(</>) 


(40) 


Indeed, gravitational theories of the Jordan-Brans-Dicke type do contain 
the possibility for a time-varying gravitational constant. However, these 
theories can always be re-expressed such that Gat is fixed and other mass 
scales in the theory become time dependent (i.e., dependent on the scalar 
field). For example, the JBD action can be written as 


S = 




T A T C' rn {^ rna tt erj 
<P 



(41) 


where ui is a number which characterizes the degree of departure from gen¬ 
eral relativity (GR is recovered as u> —> oo), A is the cosmological constant, 
and the matter action for electromagnetism and a single massive fermion 
can be written as 

C m = --^F 2 -*F 0V -mflfa. (42) 

4e z 

Written this way, if the scalar field <f> evolves, then Gn does as well. In 
another conformal frame, the JBD action can be rewritten as 


• 

(43) 

In this frame, Newton’s constant is constant, but the fermion mass (after 
T is rescaled) varies as and the cosmological constant varies as 

The physics described by either of these two actions is identical and the two 
frames can not be distinguished as the measurable dimensionless quantity 
Gm 2 oc (f >- 1 in both frames. While the fine-structure constant remains 
constant in this construction, it is straight forward to consider theories 


S = J d A x\fg 




( j ) 3 / 2 ( j ) 2 


1 

4e 2 
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where it is not. In what follows, I will restrict attention to variations in the 
fine-structure constant. 

In any unified theory in which the gauge fields have a common origin, 
variations in the fine structure constant will be accompanied by similar 
variations in the other gauge couplings 147 (see also, 148 ). In other words, 
variations of the gauge coupling at the unified scale will induce variations 
in all of the gauge couplings at the low energy scale. 

It is easy to see that the running of the strong coupling constant has 
dramatic consequences for the low energy hadronic parameters, including 
the masses of nucleons 147 . Indeed the masses are determined by the QCD 
scale, A, which is related to the ultraviolet scale, Muvi by dimensional 
transmutation: 


a s{^uv) 


9 2 s {M 2 uv ) 

47r 


47T 

6 3 ln(M2 y /A 2 )’ 


(44) 


where 63 is a usual renormalization group coefficient that depends on the 
number of massless degrees of freedom, running in the loop. Clearly, 
changes in g s will induce (exponentially) large changes in A: 


AA 

"X 


27T 


A a s {M uv ) A a s (Muv) 


9 a s (Muv) oi s (Muv) a s (Muv ) 

where for illustrative purposes we took the beta function of QCD with 
three fermions. On the other hand, the electromagnetic coupling a never 
experiences significant running from Mpy to A and thus A A/A A a/a. 
A more elaborate treatment of the renormalization group equations above 
M z 149 leads to the result that is in perfect agreement with 14 ': 


^ = 30^. 
A a 


(46) 


In addition, we expect that not only the gauge couplings will vary, but all 
Yukawa couplings are expected to vary as well. In 147 , the string motivated 
dependence was found to be 


Ah A au 

h ajj 


(47) 


where ajj is the gauge coupling at the unification scale and h is the Yukawa 
coupling at the same scale. However in theories in which the electroweak 
scale is derived by dimensional transmutation, changes in the Yukawa cou¬ 
plings (particularly the top Yukawa) leads to exponentially large changes 
in the Higgs vev. In such theories, the Higgs expectation value corresponds 
to the renormalization point and is given qualitatively by 


v ~ Mpexp(—2TTc/at) 


( 48 ) 
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where c is a constant of order 1, and a t = Thus small changes in 

ht. will induce large changes in v. For c ~ ht ~ 1, 


v ajj 


(49) 


This dependence gets translated into a variation in all low energy particle 
masses. In short, once we allow a to vary, virtually all masses and cou¬ 
plings are expected to vary as well, typically much more strongly than the 
variation induced by the Coulomb interaction alone. Unfortunately, it is 
very hard to make a quantitative prediction for Av/v simply because we do 
not know exactly how the dimensional transmutation happens in the Higgs 
sector, and the answer will depend, for example, on such things as the 
dilaton dependence of the supersymmetry breaking parameters. This un¬ 
certainty is characterized in Eq. (48) by the parameter c. For the purpose 
of the present discussion it is reasonable to assume that Av/v is compa¬ 
rable but not exactly equal to AA/A. That is, although they are both 
0(10 — 100)Aa-/a, their difference |AA/A — Av/v\ is of the same order of 
magnitude which we will take as ~ 50 A a/a. 

Much of the recent excitement over the possibility of a time variation 
in the hue structure constant stems from a series of recent observations of 
quasar absorption systems and a detailed fit of positions of the absorption 
lines for several heavy elements using the “many-multiplet” method 150,151 . 
A related though less sensitive method for testing the variability of a, is 
the alkali doublet method, which neatly describes the physics involved. 

Absorption clouds are prevalent along the lines of sight towards distant, 
high redshift quasars. As such, the quasar acts as a bright source, and the 
absorption features seen in these clouds reflect their chemical abundances. 
Consider an absorption feature in a doublet system involving for example, 
S1/2 — > P3/2 and S 1/2 —» Pi/ 2 transitions. While the overall wavelength 
position of the doublet is a measure of the redshift of the absorption cloud, 
the separation of the two lines is a measure of the fine structure constant. 
This is easily seen by recalling the energy splitting due to the spin-orbit 
coupling, 


A E ~ —• L ~ me 8 , —— ~ e 4 ~ a 2 . (50) 

Since the line splitting AA/A ~ A E/E, the relative change in the line 
splitting is directly proportional to Aa/a. The many multiplet method 
compares transitions from different multiplets and different atoms and uti¬ 
lizes the effects of relativistic corrections on the spectra. The alkali doublet 
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method 152 has been applied to quasar absorption spectra, but the sensitiv¬ 
ity of the method only limits the variation in a within an of order 10~ 5 . 
Similarly, at present, considerations based on OIII emission line systems 153 
are also only able to set limits on the variation of a at the level of 10~ 4 . 

In contrast, the many multiplet method based on the relativistic cor¬ 
rections to atomic transitions using several transition lines from sev¬ 
eral elemental species allows for sensitivities which approach the level of 
10~ 6 150 , 151 , 154 . This method compares the line shifts of elements which 
are particularly sensitive to changes in a with those that are not. At rel¬ 
atively low redshift (z < 1.8), the method relies on the comparison of Fe 
lines to Mg lines. At higher redshift, the comparison is mainly between 
Fe and Si. At all redshifts, other elemental transitions are also included in 
the analysis. Indeed, when this method is applied to a set of Keck/Hires 
data, a statistically significant trend for a variation in a was reported: 
A a/a = (—0.54 ± 0.12) x 10 -5 over a redshift range 0.5 < z < 3.0. The 
minus sign indicates a smaller value of a in the past. 

More recent observations taken at VLT/UVES using the many multiplet 
method have not been able to duplicate the previous result 154,155 . The use 
of Fe lines in 155 on a single absorber found A a/a = (0.01 ± 0.17) x 10 -5 . 
However, since the previous result relied on a statistical average of over 100 
absorbers, it is not clear that these two results are in contradiction. In 154 , 
the use of Mg and Fe lines in a set of 23 high signal-to-noise systems yielded 
the result Aa/a = (—0.06 ± 0.06) x 10~ 5 and therefore represents a more 
significant disagreement and can be used to set very stringent limits on the 
possible variation in a. 

There exist various sensitive experimental checks that constrain the vari¬ 
ation of coupling constants (see e.g., 156 ). Limits can be derived from 
cosmology (from both big bang nucleosynthesis and the microwave back¬ 
ground), the Oklo reactor, long-lived isotopes found in meteoritic samples, 
and atomic clock measurements. 

The most far-reaching limit (in time) on the variation of a comes from 
BBN. The limit is primarily due to the limit on 4 He. Changes in the fine 
structure constant affect directly the neutron-proton mass difference which 
can be expressed as Atoat ~ aah-QCD + bv, where A qcd ~ 0(100) MeV 
is the mass scale associated with strong interactions, v ~ 0(100) GeV 
determines the weak scale, and a and b are numbers which fix the final 
contribution to Atojv to be —0.8 MeV and 2.1 MeV, respectively. From the 
previous discussion on BBN, changes in a directly induce changes in Atojv, 
which affects the neutron to proton ratio. The relatively good agreement 
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between theory and observation, |AF/Y| < 3.5% allows one to set a limit 
|Aa/a| < 0.06 (A Y/Y scales with A a/a) 157 > 147 , 1 23 . Sj nce this limit is 
applied over the age of the Universe, we obtain a limit on the rate of change 
\a/a\ < 4x 10 -12 yr _1 over the last 13 Gyr. When coupled variations of the 
couplings are considered, the above bound is improved by about 2 orders 
of magnitude to Aa/a < 10~ 4 as confirmed in a numerical calculation 158 . 

One can also derive cosmological bounds based on the microwave back¬ 
ground. Changes in the fine-structure constant lead directly to changes in 
the hydrogen binding energy, Ej,. As the Universe expands, its radiation 
cools to a temperature, Td ec , at which protons and electrons can combine 
to form neutral hydrogen atoms, allowing the photons to decouple and free 
stream. Measurements of the microwave background can determine this 
temperature to reasonably high accuracy (a few percent) 9 . At decoupling 
l / _1 exp (—Eb/Tdec) ~ I- Thus, changes in a of at most a few percent can 
be tolerated over the time scale associated with decoupling (a redshift of 
z ~ 1100) 159 . 

Interesting constraints on the variation of a can be obtained from the 
Oklo phenomenon concerning the operation of a natural reactor in a rich 
uranium deposit in Gabon approximately two billion years ago. The ob¬ 
served isotopic abundance distribution at Oklo can be related to the cross 
section for neutron capture on 149 Sm 160 . This cross section depends sen¬ 
sitively on the neutron resonance energy E r for radiative capture by 149 Sm 
into an excited state of 150 Sm. The observed isotopic ratios only allow a 
small shift of |A2S r | < E r from the present value of E r = 0.0973 eV. This 
then constrains the possible variations in the energy difference between the 
excited state of 150 Sm and the ground state of 149 Sm over the last two bil¬ 
lion years. A contribution to this energy difference comes from the Coulomb 
energy Ec = (3/5)(e 2 /r 0 )Z 2 /A 4 / 3 (ro = 1.2 fm) for a nucleus with Z pro¬ 
tons and (A — Z) neutrons. This contribution clearly scales with a and is 
£c( 150 Sm) — Uc , ( 149 Sm) = 1.16(a/ao) MeV, where «o is the present value 
of a. Considering the time variation of a alone, |AU r | ~ 1.16|Aa/a| MeV 
and a limit |Aa/a| < 1CV 7 can be obtained 160 . However, if all fundamental 
couplings are allowed to vary interdependently, a much more stringent limit 
|Aa/a| < (1 — 5) x 1(V 10 may be obtained 161 . 

Bounds on the variation of the fundamental couplings can also be ob¬ 
tained from our knowledge of the lifetimes of certain long-lived nuclei. In 
particular, it is possible to use precise meteoritic data to constrain nuclear 
decay rates back to the time of solar system formation (about 4.6 Gyr ago). 
Thus, we can derive a constraint on possible variations at a redshift z ~ 0.45 
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bordering the range (z = 0.5-3.0) over which such variations are claimed to 
be observed. The pioneering study on the effect of variations of fundamental 
constants on radioactive decay rates was performed by Peebles and Dicke 
and by Dyson 162 . The /3-decay rate, A^, depends on some power n of the 
energy Qp released during the decay, A p oc Qp. A contribution to Qp again 
comes from the Coulomb energy Ec oc a. Isotopes with the lowest Qp are 
typically most sensitive to changes in a as AXp/Xp = n(AQp/Qp) is large 
for small Qp. The isotope with the smallest Qp (2.66 ± 0.02 keV) is 187 Re, 
which decays into 18 'Os. If some radioactive 18 'Re was incorporated into a 
meteorite formed in the early solar system, the present abundance of 187 Os 
in the meteorite is ( 18 'Os)o = ( 187 Os)j + ( 187 Re)o[exp(Ai 87 t 0 ) —1], where the 
subscripts “i” and “0” denote the initial and present abundances, respec¬ 
tively, Ai 87 is A p for 187 Re, and t a is the age of the meteorite. The above 
correlation between the present meteoritic abundances of 18 'Os (daughter) 
and 187 Re (parent) can be generalized to other daughter-parent pairs. All 
these correlations can be used to derive the product of the relevant decay 
rate and the meteoritic age. Using the decay rates of 238 U and 235 U from 
laboratory measurements, the correlations for the 206 Pb- 238 U and 20 'Pb- 
235 U pairs give a precise age of t a = 4.558 Gyr for angrite meteorites 163 . 
This determination of t a has the advantage that the decay rates of 238 U 
and 235 U, and hence t a , are rather insensitive to the variation of fundamen¬ 
tal couplings 162 . The above age for angrite meteorites allows for a precise 
determination of Ai 87 from the correlation for the 18 '0s- 187 Re pair in iron 
meteorites formed within 5 Myr of the angrite meteorites 164 . Comparing 
this value of Ai 87 , which covers the decay over the past 4.6 Gyr, with the 
present value from a laboratory measurement 165 limits the possible vari¬ 
ation of a to Aa/a — (8 ± 8) x 10 ' 166 . Once again, if all fundamen¬ 
tal couplings are allowed to vary interdependently, a more stringent limit 
Aa/a = (2.7 ± 2.7) x 10 -8 may be obtained. 

Finally, there are a number of present-day laboratory limits on the vari¬ 
ability of the fine-structure constant using two kinds of atomic clocks: 
one based on hyperfine transitions involving changes only in the total 
spin of electrons and the nucleus and the other based on electronic tran¬ 
sitions involving changes in the spatial wavefunction of electrons. The 
electronic transition frequency v e i depends on a relativistic correction 
Freiipt), which is a function of a and is different for different atoms. Rel¬ 
ative to v e i, the hyperfine transition frequency Vhf has an extra depen¬ 
dence on {p,nuci/H b)q?, where ii nuc i is the magnetic moment of the rel¬ 
evant nucleus and p,B is the atomic Bohr magneton. For atoms A and 
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B, Vhf,A/vhf,B CX (/^nuci,A//^nncZ,s)-^rei,A(^)/-^reZ,B(o^) find ^hf : A/^el,B (X 
(Unuci/H b)® 2 Frei^Aicz)/Frei,B(<x)■ If only the variation of a is considered, 
this can be tested by comparing Vhf,A with Vhf,B or Vhf,A with v e i t B over 
a period of time. Three recent experiments have led to marked improve¬ 
ment in the limit on the variation of a: A a/a < 6 x 10 -15 from comparing 
hyperfme transitions in 87 Rb and 133 Cs over a period of about 4 years 167 , 
A a/a < 4 x 10~ 15 from comparing an electric quadrupole transition in 
199 Hg + to the ground-state hyperfme transition in 133 Cs over a 3 year 
period 168 , and A a/a = (1.1 ± 2.3) x 10” 15 from comparing the 1S-2S tran¬ 
sition in atomic hydrogen to the hyperfme transition in 133 Cs over a 4 year 
period 169 . If both a and ^nuci/^B are allowed to vary, then constraints 
on these two distinct variations can be obtained by combining the latter 
two experiments, which give Aa/a = (—0.9 ±4.2) x 10~ 15 or d/a < 10~ 15 

yr -i 169 

A summary of the constraints on a is found in Fig. 25, taken from ref. 

170 

The result found in 154 and in the statistically dominant subsample of 
74 out of the 128 low redshift absorbers used in 151 are sensitive to the 
assumed isotopic abundance ratio of Mg. In both analyses, a solar ratio of 
24 Mg: 25 Mg: 26 Mg = 79:10:11 was adopted. However, the resulting shift in a 
is very sensitive to this ratio. Furthermore, it is commonly assumed that the 
heavy Mg isotopes are absent in low metallicity environments characteristic 
of quasar absorption systems. Indeed, had the analyses assumed only pure 
24 Mg is present in the quasar absorbers, a much more significant result 
would have been obtained. The Keck/Hires data 151 would have yielded 
A a/a = (—0.98 ±0.13) x 10 -5 for the low redshift subsample and A a/a = 
(-0.36 ± 0.06) x 10~ 5 for the VLT/UVES data 154 . 

The sensitivity to the Mg isotopic ratio has led to a new interpreta¬ 
tion of the many multiplet results 1 ' 1 . The apparent variation in a in the 
Fe-Mg systems can be explained by the early nucleosynthesis of 25 > 26 Mg. 
The heavy Mg isotopes are efficiently produced in intermediate mass stars, 
particular in stars with masses 4-6 times the mass of the sun, when He 
and H are burning in shells outside the C and O core. There may even be 
evidence for enhanced populations of intermediate mass stars at very low 
metallicity. 

Recall the dispersion seen in D/H observations in quasar absorption 
systems as seen in Fig. 18. Is there a real dispersion in D/H in these 
high redshift systems? The data may show an inverse correlation of D/H 
abundance with Si 100 ’ 102 . This may be an artifact of poorly determined 
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Figure 25. Constraints on the rate of variation (da/dt)/a as a function of the fractional 
“look-back” time (to — £)/to, where to ~ 13 Gyr is the present age of the Universe. The 
results shown are taken from data on atomic clocks 168 , considerations of the Oklo 
phenomenon 160 , meteoritic data on 187 Re decay 164 > 1 66^ anc j man y_ mu ltiplet (MM1 151 , 
MM2 154 , MM3 155 ) and alkali-doublet (AD 152 ) analyses of quasar absorption spectra. 
For convenience, the results for MM1, MM2, MM3, and AD are shown at the mean 
redshift for the data used and then converted to (to — t) /to - Note that the result for 
MM1 actually covers (to — t)/to = 0.37-0.84. Except for this result, all others are 
consistent with no time variation of a. 


Si abundances, or (as yet unknown) systematics affecting the D/H deter¬ 
mination in high-column density (damped Lyman-cc, hereafter DLA) or 
low-column density (Lyman limit systems) absorbers. On the other hand, 
if the correlation is real it would indicate that chemical evolution processes 
have occurred in these systems and that some processing of D/H must have 
occurred even at high redshift. 

It is interesting to speculate 172 that the possible high redshift destruc¬ 
tion of D/H is real and related to the chemical evolutionary history of high 
red shift systems. For example, these observations could be signatures of 
an early population of intermediate-mass stars characterized by an initial 
mass function different from that of the solar neighborhood. An example 
of such an IMF is shown in Fig. 26 173 . 

There are a number of immediate consequences of an IMF of the type 
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Figure 26. An IMF with an early enhancement of intermediate mass stars 173 . 


shown in Fig. 26. In addition to the destruction of D/H at low metallicity, 
one expects observable C and N enhancements in high redshift absorption 
systems. In addition, one also expects an enhancement of the heavy Mg 
isotopes, 25 > 26 Mg, which may account 171,173 for the apparent variation of 
the fine-structure constant in quasar absorption systems. In this sense, 
the many multiplet method can be used to trace the chemical history of 
primitive absorption clouds 1 ' 1 . This hypothesis will be tested by future 
observations and examinations of correlation among other heavy elements 
produced in intermediate mass stars 174 . 
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